arXiv:1505.01088vl [astro-ph.SR] 5 May 2015 


Astronomy & Astrophysics manuscript no. ACOROX'25734'2astroph © ESO 2017 

April 15, 2017 


Pulsations of rapidly rotating stars 
II. Realistic modelling for intermediate-mass stars 

R-M. Ouazzani^, I. W. Roxburgh^, and M-A. Dupret^ 

* Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, 
Denmark, e-mail: rhita-inaria.ouazzani@phys.au.dk 
^ Astronomy Unit, Queen Mary University of London, Mile End Road, London, El 4NS, UK 
^ Institut d’Astrophysique et de Geophysique de I’Universite de Liege, Allee du 6 Aout 17, 4000 Liege, Belgium 

Received January 21st / Accepted May 4th 


ABSTRACT 

Context. Very high precision seismic space missions such as CoRoT and Kepler provide the means for testing the modelling of 
transport processes in stellar interiors. For some stars, such as S Scuti y Doradus and Be stars, for instance, the observed pulsation 
spectra are modihed by rotation to such an extent that it prevents any fruitful interpretation. 

Aims. Our aim is to characterise acoustic pulsation spectra of realistic stellar models in order to be able to interpret asteroseismic data 
from such stars. 

Methods. The 2-dimensional oscillation code ACOR, which treats rotation in a non-perturbative manner, is used to study pulsation 
spectra of highly distorted evolved models of stars. 2D models of stars are obtained by a self-consistent method which distorts 
spherically averaged stellar models a posteriori, at any stage of evolution, and for any type of rotation law. 

Results. Four types of modes are calculated in a very dense frequency spectrum, among which are island modes. The regularity of 
the island modes spectrum is confirmed and yields a new set of quantum numbers, with which an echelle diagram can be built. Mixed 
gravito-acoustic modes are calculated in rapidly rotating models for the first time. 
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1. Introduction 


Rapid stellar rotation introduces a number of phenomena which 
considerably complicate the modelling of stars and the study of 
their pulsation modes. On the one hand, centrifugal accelera¬ 
tion reduces local gravity, mimicking a lower mass. On the other 
hand, rotation induces meridional circu lation (|Eddington|1925| l 
and shear and baroclinic instabilities ( [Mathis & Zahn||2004| l, 
which contribute to the mixing of chemical elements. Although 
crucial for modelling stellar structure and evolution, the problem 
of transport of angular momentum inside stars has not yet been 
fully understood. 

As shown with great success for evolved low mass stars, 
asteroseismology is a very efficient tool to constrain the inter¬ 
nal structure of stars. Unfortunately, for intermediate mass stars, 
because they do not undergo magnetic braking, rapid rotation 
comes into play and makes the interpretation of their seismic 
spectra much more delicate. East rotation is known to have a 
strong impact on stellar pulsations: the centrifugal force distorts 
the pulsations cavity and the Coriolis force modifies the dynamic 
of modes. So far the impact of fast rotation on observational fre¬ 
quency spectra is poorly understood, and some questions still 
need to be answered. What are the modes corresponding to each 
observed frequency? Which cavity inside the star do they probe? 
How can the physics inside this cavity be extracted? All these 
questions call for realistic modelling. 

Eor slow rotators, the effect of rotation are taken into account 
as a small perturbation of pulsations (see for instance jLedouxj 
|1951| for first-order effects, the second-order effects have first 
been explored by Vorontsov|1981| 1983 see also Dziembowski 
|& Goodel 19921 and we refer to jSoufi et al.|1998| for third-order 


effects). But these methods, although elegant and simple to use, 
have their limitations (Reese et al. 2006 [ Suarez et al. 20101; they 
fail to adequately model the effects of rapid rotation in stars that 
show very high surface velocities, such as 6 Scuti and Be stars, 
nor the impact of the Coriolis force on gravity modes in stars 
whose surface rotates slowly, but in which the pulsation periods 
are of the same order as their rotation period (Pror ~ Pose), such 
as SPB or y Doradus stars 


This is important for many stars in the CoRoT and Kepler 
fields of observation. Hence, if one aims at correctly extract¬ 
ing the internal structure from seismic observations, they need 
to correctly describe the effects of rotation on pulsations. In that 
context, the aim of this article is to investigate synthetic seismic 
spectra of such stars for any type of rotation profile Q(r, 6), dif¬ 
ferential both in radius r and co-latitude 6, and to propose new 
seismic diagnostics in order to interpret their seismic observa¬ 
tions. 


In a previous paper Roxburgh P006|l presented self- 


consistent 2-dimensional models of main-sequence stars where 
the angular velocity is differential both in depth and in latitude, 
and for different stages of evolution. The adiabatic oscillations of 
stars are only governed by the acoustic structure of the star, that 
is the hydrostatic structure; density, pressure, gravity and adia¬ 
batic exponent. The pulsations do not depend on whether the star 
is in thermal equilibrium. The self-consistent method allows one 
to determine such an acoustic structure. 


In the first paper ( jOuazzani et al.|[20I2| hereafter Paper 
I), we presented the ACOR pulsation code -Adiabatic Code 
of Oscillation including Rotation-, and its results for a 2- 
dimensional polytropic model. Here, we report on pulsation 
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modelling of an evolved distorted model of an intermediate mass 
star, typical of a 5 Scuti star. Section introduces the method 
applied to obtain distorted models of stars following [Roxbur^ 
( |2006| l. Section is dedicated to the pulsation modelling with 
an emphasis on how the ACOR code was improved in order to 
compute the pulsations of realistic distorted stellar models. In 
Section we present an application to a typical model of in¬ 
termediate mass star. Finally, Sect, [^reports on regularities re¬ 
vealed in the theoretical pulsation spectrum of the model, and 
conclusions are given in Sect. 


2. Self-consistently distorted model 

Stellar evolution in two dimensions has yet to be modelled. 
Progress is being made in this field for instance by the ESTER 
project (for Evolution STEllaire en Rotation, Rieutord [2007 


Espinosa Lara & Rieutord] 2013 


. The ultimate goal of the 
ESTER project is to solve the hydrodynamics equations in 
two dimensions along with the chemical evolution. Currently 
ESTER can only model the steady state with regards to evolu¬ 
tion, the central chemical composition being specified ad hoc. 

An alternative to this approach is to calculate spheri¬ 
cally averaged stellar evolution models, then distort the mod¬ 
els a posteriori using a self-consistent method. The first self- 
consistent model was developed for rapidly rotating polytropes 
and white dwarfs by [James ( 1964|). See also Hurley & Roberts 


(|1964|l, closely followed by a two-zone model by [Monaghan 


& Roxburgh (19651 and two-zone models for rapidly rotat¬ 
ing main-sequence stars by Roxburgh et al. ( 1965[ ). In paral¬ 
lel, [Ostriker & Mark| ( [I968| l had developed and used such an 
approach to compute uniformly and differentially rotating poly¬ 
tropic stellar models. The method has then been improved by 
including more realistic input physics Packson|1970| l, but only 
applied to massive stars, as convergence problems arose when 
dealing with low mass stars, due to the high concentration of the 
central mass ( [Clement|1978| l. [Jackson et al.[ ( |2004) l; |MacGregor| 
et al. (2007 1 solved this problem and successfully applied the 
method to intermediate mass stars. In their models, the only re¬ 
maining limitation is the chemical homogeneity and conserva¬ 
tive laws of rotation. 

We present here such an approach based on a method es¬ 
tablished by |Roxburgh| ( [2006] ) which allows the computation of 
non-homogeneous, fully differentially rotating stellar models. It 
consists of building the two-dimensional acoustic structure of 
a rotating star in hydrostatic equilibrium starting with spherical 
profiles of the structural quantities using an iterative scheme. 


2.1. Acoustic structure 


We recall that the adiabatic pulsations of stars are solely gov¬ 
erned by their acoustic (i.e. hydrostatic) structure, that is the in¬ 
ternal density p, pressure p and gravity Vd), together with the 
adiabatic exponent Ei. As shown in Roxburgh[ ( 2006[ l, for rotat¬ 
ing stars the 2-dimensional hydrostatic structure p(r, 6), P{r, 6) 
and gravitational potential <I>(r, ff) are determined solely by the 
angular velocity profile Q(r, 6), the density p„,(r) = p(r, 6„,) and 
the surface pressure Ps along any one angle 0^. The adiabatic 
exponent ri(r, 0) is determined from p(r,6), P(r,6) and the 2- 
dimensional hydrogen profile X{r, 0). If we assume a particu¬ 
lar model of chemical elements transport (which is the case in 
this work: X constant on surfaces of constant entropy) then it 
is sufficient to have X„,{r), along 6„,. The method then consists 
of generating acoustic models of rotating stars with any Q.(r, 0) 
by considering that the radial density and hydrogen abundance 


along the angle 0„, take the values from a spherically averaged 
model, yielded by 1-dimensional evolution modelling. The ref¬ 
erence angle 6m is taken here as the root of the second order 
Legendre polynomial: 6m = cos“'(l/ V3). 

The scheme is the following: 

1. Take as an initial guess p(r, 0) = Pm{r) 

2. Solve the Poisson equation: 

V^O = 47rGp (1) 


for <l>(r, 0) given p(r, 0), together with the proper boundary con¬ 
ditions. This is done by expanding <l>(r, 0) and p(r, 0) into series 
of Legendre polynomials. The density decomposition is found 
by imposing that p(r, 0) = p{ri, 6j) at each grid point (r,-, dj). 
Given that, the coefficients of the gravitational potential expan¬ 
sion are determined from an integral representation of the solu¬ 
tion of the Poisson equation and the upper boundary condition 
projected onto the Legendre polynomials (for a detailed solution, 
see |Roxburgh|2004| ). 

3. Solve the curl of the hydrostatic equilibrium: 


curl 


VO-Q^ 


= 0 


( 2 ) 


for p(r,6) given 0(r, 0) from step 2. m - (r sin 0, r cos 0,0) is 
the vectorial distance from the rotation axis. This is done in 
two steps: first, we seek for the characteristics curves r(0) along 
which the partial differential equation is an ordinary differential 
equation (ODE), and then solve E q. ([2 1 as an ODE along these 
curves (for a detailed solution, see Roxburgh[2006 1. 

4. Then this density profile is injected at step 2 and steps 2 
and 3 are iterated until convergence of the equatorial and the po¬ 
lar radii. 

5. The characteristic curves turn out to be the isobars. Once 
the characteristic curves are determined, all we need is the value 
of P along a radius at any angle 0, P(r, 6) is then constant along 
the characteristic curves. 

6. To determine the 2-dimensional profile of the adiabatic 
exponent given P{r, 6) and p(r, 0), we assume that Xir, 0) is con¬ 
stant on isobars (justification for this assumption follows in Sect. 
\2.2) . We then solve the equation of state for Ti = Yx{P,p,Xm), 
Xm being the value of the hydrogen mass fraction at the reference 
angle 0,„. This gives Ei along the isobars. For further details we 
refer to Roxburgh (20061. 


2.2. Discussion 

For the algorithm to converge, several approximations must be 
assumed, the first of which concerns the structure of surface lay¬ 
ers. In stars that are differentially rotating, the fluid is baroclinic. 
This means that surfaces of constant pressure do not correspond 
to surfaces of constant density. Hence the solution of the algo¬ 
rithm presented in the previous section may have p = 0 on sur¬ 
faces of P > 0. For the solution to be well behaved at the surface, 
we need curl(Q^^) = 0 in the surface layers. This is ensured by 
taking Q. constant in the outer layers of the star. 

Moreover, in this region of low density and pressure, the 
meridional circulation driven by the rotation can become very 
large, and induce small scale turbulence. We assume here that 
the time scale of the circulation is very small in these layers, 
so that transport of angular momentum is rapid, and the system 
reaches a steady state (Q constant near the surface). 

Finally, in order to retrieve the adiabatic exponent Fi, we 
have assumed that the hydrogen mass fraction X is constant 
on isobars. However, in the fully ionised interior of most stars 
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Fig. 1. Coordinate system used in ACOR. Vq corresponds to a 
region inside the convective core, V i extends from a boundary 
inside the convective core to the stellar surface, and V 2 encom¬ 
passes the star. 


where, due to nuclear reactions, X differs from its initial value 
Xq, the values of F 1 are only weakly dependent on X. On the 
other hand in the outer layers, where F 1 can vary considerably, 
in absence of diffusion X is unchanged from its initial value Xq. 
Therefore X constant on isobars seems to be a reasonable as¬ 
sumption here. 


3. Pulsations modelling 

As mentioned above, for the computation of pulsations we only 
require the hydrostatic structure of the star and not the thermal 
structure. We compute the oscillation modes as the adiabatic re¬ 
sponse of the structure to small perturbations - i.e. of the den¬ 
sity, pressure, gravitational potential and velocity field - using 
the Eulerian formalism. 

The code ACOR has been developed for this purpose and 
is presented in Paper 1. Here we report on new improvements 
which were necessary to implement in order to compute the pul¬ 
sations of non-barotropic models of stars, which present steeper 
gradients around the poles than near the equator. 

3.1. A new coordinate system 

The method reported in Sect. yields a 2-dimensional stellar 
model described on the characteristics curves (isobars). In order 
to avoid potential uncertainties caused by interpolation, the 
idea is to keep the structural quantities expressed on the isobars 
as much as possible when computing the pulsations of such 
models. Nevertheless, to simplify boundary conditions -at the 
centre, where the eigenfunctions should behave spherically, and 
at inhnity, where the gravitational potential should vanish- a 
new hybrid coordinate system is proposed. The new radial coor¬ 
dinate ^ is expressed as a bijective, monotonic and continuously 
differentiable function of the radius r. From spherical coordi¬ 
nates at the centre, the system progressively matches isobars, 
and from the distorted stellar surface becomes spherical again 
outside the star. Therefore, we adopted a multidomain approach, 
which consists in dividing the physical space into 3 domains 
whose boundaries correspond to the models discontinuities: 
one central domain Vq whose upper boundary is included in 
the convective core, one domain Vi which encloses the stellar 


surface, and one external domain V 2 , bounded by the stellar 
surface and an external sphere. 


Domain Vq from { — Oto ^ - Xc,eq 

Here Xc,eq is the ratio of the radius at the edge of the convective 
core over the total radius of the star, both at equator (see Fig. 

the inner domain, encircled by the red boundary). In Paper 
1, where polytropic structural models were treated, the bijective 
transformation continuously goes from spherical at the centre to 
matching the stellar surface without intermediately adopting the 
isobars. In the case of a non-barotropic model of star, isobars are 
closer to each other at the poles than at the equator. To ensure 
continuity of d^r over the boundary ^ = Xc,eq, d^r should de¬ 
pend on the latitude when ( goes to Xc.eq. Hence, the following 
transformation has been developed: 


V^c =^/Xc,eq € [0;1] ,V0 e [0;;r] , 

x(^,e) = 


(3) 


^c,eq 


543 - 34^ (xM , 
2 Uc,eq 


^c,eq 


d( 


where x is the fractional radius (r/R), and dxdd^ is the 
derivative of x with respect to 4" at f = Xc,eq. The first term 
in brackets enforces a spherical behavior for x{^, 0) near the 
center, the second term allows to match the isobar Xc(0) located 
at = 1 and the last term insures that the derivative of x(^, 6) at 
= 1 matches d^xdO). We chose the simplest functions which 
fulhlled the constrains. Odd polynomials of d have been chosen 
because they go faster toward zero at the centre. In practice, 
Xc,eq can take any value smaller than the actual fractional radius 
of the fj. gradient location, both for low mass and massive stars. 


Domain Vi from ^ = Xc.eq to ^ = 1 

In this domain, constant ^ correspond to the isobars, with 
( = x(6„,). x(d, 0) is determined by the characteristics equation. 


Domain V 2 from ^ = 1 to f = 2 

In 1^2 the coordinate system goes gradually from the distorted 
stellar surface to a sphere. A new transformation is proposed, 
which ensure the continuity of r and d^r at the surface of the 
star ^ = 1. The same mathematical form as for Eq.([^ applies 
here, except that Eq.Q is not taken in zero, therefore, there is 
no need for odd polynomials. The simplest form which fulfills 
the constrains is: 


6 [1;2] ,V6» e [0;7r] , (4) 

xid 0)^i^ + (2^3 - 9^2 + 12^ - 4) (x,(0) - 1) 

+ (^3_5^2^8^-4) J^(0)-lj 

In the two domains Vq and V 2 where the coordinate system 
does not rely on isobars, the structural quantities are interpolated 
on the new grid (f, 0) by the help of cubic polynomial interpola¬ 
tion over two grid points. 

3.2. The pulsations computations 

It consists of solving the hydrodynamics equations perturbed by 
Eulerian fluctuations, by direct integration of the 2-dimensional 
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problem. The numerical method is based on a spectral multi- 
domain method which expands the angular dependence of eigen¬ 
functions into spherical harmonics series, and whose radial treat¬ 
ment is particularly well adapted to the behaviour of equilibrium 
quantities in evolved models (at the interface of convective and 
radiative regions, and at the stellar surface). The radial differ¬ 
entiation is made by means of a sophisticated finite difference 
method which is accurate up to the fifth order in terms of the 
radial resolution (developed by Scuflaire et al.|2008) l. 

This code has been validated by comparison with the results 
of Reese et al. ( |2006| ) for polytropic models. The agreement be¬ 
tween the two codes is found excellent. All details concerning 
the ACOR code are given in Paper 1. 


4. An example, M - 2Mq, Xc - 0.35, O = Q.{r) 

4.1. The stellar model 

As an illustrative example we take Q = Tls(^(x) as a function 
only of radius x = r/Ro, where a>{x), shown in Fig. decreases 
from a value of 3 in a central core to 1 in the outer envelope, 
and has continuous derivatives. As mentioned in Sect. |2.2| with 
Q constant in the outer layers the solution is well behaved at the 
surface. For the spherically averaged model we took a star of 2 
Mq with an initial composition X = 0.72, Z = 0.02, evolved to 
the stage where = 0.35. The average centrifugal force 2Q^r/3 
was added to the hydrostatic equation, the remaining equations 
being those of spherical star. The dimensionless angular velocity 
a>{x) was taken as fixed throughout the evolution and Qj deter¬ 
mined by requiring conservation of angular momentum for the 
star as a whole. The initial model had radius R ^ 2.5 Rq and 
Q.S = 4.3 X lO^^rad.s^^ corresponding to an equatorial veloc¬ 
ity of ^ 80km.s '. The star was evolved to a central hydrogen 
abundance X^. = 0.35, at which stage flj = 1.296 x lO^^rad.s^'. 

The models were produced using the STAROX evolution 
code (cf. Roxburgh|2008[) wi th the OPAL2001 equation of state 
( Rogers & Nayfonov||2002|), OPAL GN93 [Iglesias & Rogers] 
( 1996|l and" Alexander & Ferguson ( |1994| l opacities, NACRE 
( Angulo et al.|T999| l nuclear reaction rates, and a radial mesh 
in mass with N = 2000 points. The spherically averaged model 
was taken to represent the 2-dimensional model along t he an gle 
6m - cos“'(l/ V3), and the procedure described in Sect. 


2.1 


_ fol¬ 
lowed to produce the 2-d model. The angular mesh had Ny = 240 
points and the mesh in domain V 2 for solving Poissons equation 
had Ni = 500 points. The flatness of the model 1 - Rpoi/Req - 
0.25, the equatorial rotational velocity Veg = 282.3km.s“*, and 
the ratio of centrifugal force to gravity at the surface equator 
Q2 r3 /GM = 0.65. 


4.2. The pulsation modes 

Given the rotation velocity (Q = 0.8 IQ^^ at the surface, where 
Qi = ^GM/Rlg) a convergence study found that 40 spherical 
harmonics were necessary in the spectral expansion of eigen- 
modes. In the present article we confine our analysis to axisym- 
metric modes, although there is no restriction to such modes in 
the ACOR code. Odd and even modes with m = 0 are presented 
here. The parity of a mode p is given by p = (^ - 1 - \m\) mod 2. 
Therefore modes with p - Q are called even modes, and modes 
with p - \ odd modes. For axisymmetric modes, odd modes 
(resp. even) are expanded on spherical harmonics series with 
only odd (resp. even) values of angular degree {. 



0 0.2 0.4 0.6 0.8 1 

C 

Fig. 2. Radial profile of the rotation angular velocity scaled by 
the equatorial angular velocity at the surface. Q is constant on 
isobars but differential in the normal direction. Therefore, the 
model is not barotropic. 



Fig. 3. Simulated spectrum of axisymmetric modes (only) for 
the model presented Sect. it contains the eigenfrequencies 
of odd and even Island modes (in blue and red resp.), as well 
as modes trapped only in the p gradient. Chaotic modes and 
Whisperring gallery modes in grey. 


In this study, we have computed axisymmetric eigenmodes 
with frequencies ranging from 250 pHz to 750 pHz (10-30 G^.), 
by scanning the frequency interval we made sure that all the ax¬ 
isymmetric modes with contributions up to 7” = 80 are found. 
The resulting pulsation spectrum presents around 800 modes, 
which corresponds to one mode every 0.63 pHz in average. It is 
shown in Fig. where we have simulated each frequency peak 
by a Lorentzian of width 0.05 pHz and height 1. The quantity on 
the y axis of that graph has no significance. The different ampli¬ 
tudes are to be attributed to the fact that when two frequencies 
are closer than 0.05 pHz, their amplitudes add up. 

Fast rotation drastically modifies the geometry of the pulsa¬ 
tion modes; their geometric structure determines their visibility, 
as well as the region they probe. Therefore, a correct character¬ 
isation of the geometry of different modes is an important issue 
for the asteroseismology of fast rotators. The geometry of each 
of the 800 modes have been visualised (based on their kinetic en¬ 
ergy in a meridional plane) allowing a careful characterisation. 
Depending on the region of the frequency spectrum explored by 
our computations, four geometric classes are found. 


4 



















































































R-M. Ouazzani et al.: Pulsations of rapidly rotating stars 



Fig. 4. Kinetic energy p\vpuis?' in a meridional plane, of an ^ = 1 
axisymmetric island mode. 



Fig. 5. Kinetic energy divided by the square root of the den¬ 
sity and multiplied by the square of the distance to the axis in 
a meridional plane, for an axisymmetric mixed island mode of 
1 . 


At low frequency, we mostly find modes trapped in the re¬ 
gion of sharp variation of the mean molecular weight, charac¬ 
terised by a sharp feature in the Brunt-Vaisala frequency. These 
modes are counterparts of high order g-modes in the rapidly ro¬ 
tating case. As in the non-rotating case, their wave number is 
proportional to yjtit + 1). Therefore, since we include compo¬ 
nents with i up to 80 here, g-modes with high ( components are 
very numerous, even in the high frequency range. Their ampli¬ 
tude is confined around the central region, they have negligible 
amplitude in outer layers, and we expect them to undergo radia¬ 
tive damping, therefore, they should not be detected observation- 
ally. 

At higher frequency, we find whispering gallery modes, 
chaotic modes and island modes such as those computed 
in homogeneous models of stars by Lignieres et al.| ( |20061 l. 
Whispering gallery modes are the counterpart of high degree 
acoustic modes in the rapidly rotating case. They probe the outer 
layers of the star, but as they show very low visibility factors, 
they might not be detected. 

Chaotic modes have no counterparts in the non-rotating case. 
They have been computed only in very high rotating models 
of stars ( Lignieres & Georgeot|2009 l. Their spatial distribution 
does not show a simple symmetry. They have significant am¬ 
plitude in the whole of the stellar interior. The lack of symmetry 
should result in a low cancellation factor. Therefore, these modes 
are expected to be detected observationally. 

Whispering gallery and chaotic modes correspond, respec¬ 
tively, to high and intermediate values of I - \m\ at zero rotation 
( Lignieres et al.|2006| l. 

Island modes are counterparts of low degree acoustic modes 
in the rapidly rotating case. They probe outer layers of the star, 
and present good geometric visibility factors. Therefore they 
should be easily detected observationally. The geometry of such 
a mode is presented Figj^with the distribution of kinetic energy 


of the mode in a meridional plane. As mentioned in Reese et al. 
(2008|l if their non-rotating counterpart is characterised by the 


known quantum numbers the geometry of these Island 

modes follows a different distribution characterised by a new set 
of quantum numbers {h - 2n + p,{ - m), where p indi¬ 

cates the mode’s parity. The two new quantum numbers € and h 
correspond respectively to the number of nodes along and per¬ 
pendicular to the ray path (Fig0, see |Reese et al.| ( |^08[|2009bl l. 

Some island modes in the explored frequency range are 
found to be mixed, i.e. they behave as pressure modes in the 


envelope and gravity modes near the core. In Figj^we give the 
geometry distribution of a mixed mode in a meridional plane, il¬ 
lustrated by the kinetic energy divided by the square root of the 
density and multiplied by the square of the distance to the axis. 
Mixed non-radial modes have been theoretically known since 
Osaki ( 1975)1 discovered them in a non-rotating model of a 10 


Mq star on the main sequence. For the first time, we show that 
such gravito-acoustic mixed modes are present in the theoreti¬ 
cal spectrum of rapidly rotating stellar models too. Indeed, as 
a star evolves, the Brunt-Vaisala frequency, i.e. the maximum 
frequency of trapped g-modes, increases owing to the increas¬ 
ing core density. Therefore, the frequencies of g modes become 
comparable to those of p modes, and some of them adopt a 
mixed character. These mixed modes are very useful because, 
unlike pure g modes, they have large enough surface amplitudes 
to be detected and yet are sensitive to the structure of the core. 
In the space photometry era, these mixed modes have been in¬ 
tensively observed in evolved solar type stars, first with CoRoT 
(|Deheuvels et al.||2010|l and with Kepler (e.g. |Campante et al. 
2011[|Mathur et al.||2011| l. These modes enabled one to probe 


the structure of the core of subgia nts (|Deheuvels & Michel 201 l| l 
and red giants ( Beck et al.|201 l[|Mosser et al.|201 1 1, thus mak¬ 
ing it possible to discriminate between evolutionary scenarios 


(shell hydrogen or core helium burning, see for instance Mosser 
et al.|201^. The detection of mixed modes that are split by ro¬ 


tation also allow to probe the rotation rate of stars even in their 
deepest interior ( [Beck et ar]|2012[ Deheuvels et al.||2012| l. The 
situation of fast rotating intermediate mass stars is much less 
favourable, not only because they are not solar-like pulsators, 
but also because the impact of rapid rotation on their spectra was 
poorly modelled. Once stellar and pulsations models are adapted 
to such impact, the use of mixed modes as probe for rapidly ro¬ 
tating stars interior seems now reachable in the near future. 


5. Island modes regularity 

Low £ island modes are expected to be the most visible modes 
in the seismic spectra of rapidly rotating stars since they are 
the rotating counterparts to modes with low € values. Their fre¬ 
quency organisation has been studied in Lignieres et al. (|2006|l 
in non-distoited polytropic models (of index n=3), [Lignieres & 
Georgeot (|2008| |2009|l u sing ray dynamics in polytropic mod- 
els, |Reese et al.| (|20081l in fully distorted polytropic models. 


[Reese et al.| (|2069a[) for homogeneous models distorted using 


5 
































































R-M. Ouazzani et al.: Pulsations of rapidly rotating stars 



frequency modulo A~ (f2|^) 


Fig. 6. Echelle diagram for axisymmetric Island modes of dif¬ 
ferent degree h from left to right, { - 2,0 and 1, for odd (blue) 
and even modes (red). 


the SCF method ( [MacGregor et al.|2007| l, and jPasek et al.| ( |2012j l 
who have developed an semi-analytic formula for their spectrum 
under the Cowling approximation, and neglecting the Coriolis 
force. 

The hrst difficulty with studying the spectrum of island 
modes is to be able to identify them among all the other type 
of modes present in the spectra of rapidly rotating stars (chaotic 
and whispering gallery modes). This was done by hnding a few 
island modes in the high frequency regime, identifying their spa¬ 
tial geometry (i.e. h, {, m is fixed) determining their frequency 
spacings and then checking visually the modes which frequen¬ 
cies are around the expected one. 

5.1. Asymptotic frequencies 

Reese et ah] ( |2008| ) suggested a new asymptotic formula for is¬ 
land modes: 

‘^h,e,m - riAn + (A.( + A.,n + a, (5) 

which they have tested for the SCF homogeneous models. Based 
on our synthetic spectrum of axisymmetric modes we were able 
to determine the parameters A,-;, and a of this asymptotic 
formula based on a fit on our frequency spectrum. This ht has 
been carried out for the fully rotating 2d model, as well as for 
the spherical non-rotating model corresponding to the structural 
quantities of the 2d model at the latitude 0„ - cos“*(l/ V3). The 
results are shown in Table ([T]l. 


Table 1. Values of the parameters of Fq. (|^ for the non-rotating 
model (first row), and for the 2D model (second row). 


Q. 

An 

A,' 

a 

at 

at 


at 

0.000 

1.184 

0.523 

1.058 

0.807 

0.795 

0.459 

1.015 


In order to highlight this regularity, one can also plot, what 
we call an echelle diagram which is the spectrum folded about 
the spacing A;;. The result is shown in Fig. which gives the 
echelle diagram for axisymmetric island modes of different par¬ 
ity and for (-0,1, and 2. For a given parity, the modes frequen¬ 
cies line up along ridges of given ( values. 


This feature is of course not surprising as island modes are 
acoustic modes, and once the shape of their cavity is known, 
their spectrum is determined by the stratification of the fluid in 
the cavity. 

As pointed out in non-rotating models of evolved low mass 
stars, when acoustic modes acquire a gravity contribution dur¬ 
ing an avoided crossing for instance, their frequencies have the 
tendency to deviate from the regular pattern defined by the large 
spacing. That is what is causing the deviation at low frequency in 
Fig.§ where the corresponding eigenmodes clearly show trap¬ 
ping in the p gradient region. For instance, the lowest frequency 
mode on the odd ( - \ ridge is the mixed mode whose kinetic 
energy distribution is given in Fig|^ 


5.2. Theoreticsi Large Separation 


In previous works of full modelling ( [Reese etar]|2008| [2009b[ l 
the hypothesis that the pseudo-large separation Apt is a measure 
of the travel time in the island modes cavity has been mentioned 
but never actually tested. The difficulty here is to find the ray 
path for island modes geometry. The approach adopted consists 
in finding the underlying ray path for our model based on the 
topology of the 2-dimensional eigenfunctions. For the sake of 
simplicity, we assume that the underlying ray path corresponds 
to the local maxima of the pressure perturbation for w(-0 
island mode. This yields the p{^) curve plotted in red in Fig. 
where p = cos 9 and f the pseudo-radial coordinate. The theo¬ 
retical large separation is then given by: 



where Cs is the sound speed along the ray, and ds corresponds 
to the arclength along the ray. Such a method allowed us to 
determine a theoretical large frequency separation of Af™ - 
0.829 Qi:, which is compatible with the large separation com¬ 
puted from the theoretical seismic spectrum by averaging the 
individual values A^ = (0.795 + 0.012)0^;, but not within the 
dispersion. 

This is a rough estimate of the inverse of the travel time 
on the island modes ray paths, but this confirms that this new 
spacing A,~, is an equivalent of the large separation Av of acous¬ 
tic modes without rotation. Once island modes are identified 
in the spectrum. A;; allows one to constrain the sound speed 
profile in the outer radiative envelope of these rapidly rotating 
intermediate-mass stars. The other way around, in theoretical 
studies, it can also be of great help in order to track island modes 
in synthetic spectra by pointing toward the right pattern. 

Most of the uncertainty of this method lies in the determi¬ 
nation of the ray path from the eigenfunctions of the modes. Of 
course the most accurate way to compute the ray path is actu¬ 
ally by studying the dynamics of acoustics rays (as is done in 
Lignieres & Georgeot 2009 [ for polytropes), but so far it has 


never been done in a realistic model, and without making certain 
assumptions such as that the Coriolis force doesn’t affect too 
much acoustic-ray dynamics ( Reese et al.|2006 have shown that 
its effects even on relatively low radial orders are weak), and the 
Cowling approximation. The presence of mixed modes in our 
spectrum of a realistic mid main sequence model requires one to 
drop these assumptions. We leave this study for future work. 


6. Conclusions 

A new modelling chain which computes non-radial adiabatic 
pulsations in rapidly rotating distorted models of stars has been 


6 





























R-M. Ouazzani et al.: Pulsations of rapidly rotating stars 


1 



0 1 


C 

Fig. 7. Kinetic energy distribution of an ^ = 0 axisymmetric 
island mode in the {(, fi) plan, where ^ is the pseudo-radial co¬ 
ordinate, and fj. is related to the colatitude by: /t = cos 6. The 
red curve connects the local maxima of the kinetic energy, that 
defines the ray path of island modes in the explored model. 


presented. The first step consists in building a distorted stellar 
model starting with a spherically averaged one for arbitrary ro¬ 
tation laws ( |Roxburgh]|2006[ ) . The second step is the computa¬ 
tion of its stellar pulsations by fully acounting for the impact of 
centrifugal distortion and for the effect of the Coriolis force on 
mode dynamics (using the 2D non-perturbative oscillation code 
ACOR, see Paper I). In rotationally distorted stars, the struc¬ 
tural gradients are steeper near the poles than near the equa¬ 
tor. To handle that, a new coordinate system was developed that 
greatly simplifies the computations, and disminishes numerical 
sources of uncertainties. That allowed to compute for the first 
time the synthetic spectrum of a realistic distorted model of a 
non-homogeneous star, evolved on its main sequence. 

This approach has been applied to a 2 Mq star evolved mid¬ 
way on its main sequence {Xc - 0.35). The simulations yield 
a wealth of modes with a spectrum denser than its non-rotating 
counter-part. Four types of modes have been computed, among 
which island modes that are expected to be the most visible ones. 
The regularity of island modes is confirmed in the general case 
of non-barotropic stratified stellar models. In particular the spec¬ 
trum is organised around the large spacing A^, and it is possible 
to build echelle diagrams for island modes. Once the shape of the 
cavity is correctly defined, it is possible to retrieve a theoretical 
value of this large spacing as the inverse of the travel time along 
the modes path which is compatible with the value computed 
from the spectrum. 

Moreover, we have shown that, just like their non-rotating 
counterparts, the island modes can be involved in avoided cross¬ 
ings with modes trapped near the core. Such mixed modes (in 
the sense of gravito-acoustic modes) deviate from the regular 
pattern defined by the large spacing A;;. 

Finally, several seismic diagnostics based on these island 
modes arise. The large spacing of the island-mode spectrum is 
related to the stratification of the star in the outer radiative en¬ 
velope. The deviation from the regular spacing due to avoided 
crossings should allow to constrain the evolutionnary status of 
the star. 
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